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We study effects of the internal coherent (Rabi) coupling in vortex lattices in two-component 
BECs under rotation. We find how the vortex lattices without the Rabi coupling known before 
are connected to the Abrikosov lattice of integer vortices with increasing the Rabi coupling. We 
find that 1) for small Rabi couplings, fractional vortices in triangular or square lattice for small or 
large inter-component coupling constitute hexamers or tetramers, namely multi-dimer bound states 
made of six or four vortices, respectively, 2) these bound states are broken into a set of dimers at 
intermediate Rabi couplings, and 3) vortices change their partners in various ways depending on the 
inter-component coupling to organize themselves for constituting the Abrikosov lattice of integer 
vortices at strong Rabi couplings. 



Multi-component condensations are one of growing 
topics in condensed matter physics such as exotic super- 
conductors, superfluid 3 He, multi-component or spinor 
Bose-Einstein condensates (BECs) of ultra-cold atomic 
gases, excit on-polar iton condensates, nonlinear optics, 
and nonlinear sciences. They also appear in high en- 
ergy physics and astrophysics such as hadronic matter 
composed of neutron and proton Cooper pairs relevant 
for cores of neutron stars, and quark matter composed 
of di-quark condensates consisting of quark Cooper pairs 
which might be present at more high density. One of 
new common features of these systems is the existence of 
exotic vortices created by rotating superfluids or BECs 
or by applying magnetic field on superconductors. There 
exist different vortices winding around different conden- 
sates in general. Their quantized circulations for super- 
fluids or BECs and fluxes for superconductors are not 
integer valued anymore but are rational or fractional in 
general, as found in various systems: superfluid 3 He [TJH], 
p-wave superconductors [U[3H5], multi-gap superconduc- 
tors [6H9], spinor BECs [El HQ, multi-component BECs 
[HHH], exciton-polariton condensates [22 , 23J, nonlinear 
optics [24J, and color superconductors as quark matter 

Among various condensed matter systems admitting 
vortices, BECs of ultra-cold atomic gases provide a par- 
ticularly ideal system for examining properties of vortices 
both theoretically and experimentally [26J. Theoretically, 
BECs can be quantitatively well described in the mean- 
field theory, i.e., the Gross-Pitaevskii (GP) equation. On 
the other hand, BECs are quite flexible and controllable 
systems experimentally, because the atomic interaction is 
tunable through a Feshbach resonance [27J and the con- 
densates can be visualized directly by optical techniques. 
Two-component BECs have been realized by using the 
mixture of atoms with two hyperfine states of 87 Rb [28J 
or the mixture of two different species of atoms [29H3T] , 



One of fascinating features of fractional vortices is 
the possibility of various structures of vortex lattices, 
as found in a vortex phase diagram in two-component 
BECs [T3HT7] . With increasing the inter-component cou- 
pling, one obtains from an Abrikosov's triangular lat- 
tice of fractional vortices to a square lattice of fractional 
vortices [T3HT5] , and a vortex sheet [19]. One may ex- 
pect a similar phase diagram in exotic superconductors 
if the rotation speed is replaced with an applied mag- 
netic field. However, for multi-gap superconductors, the 
existence of a Josephson coupling between different con- 
densates is inevitable. While this term provides a gap 
to the Leggett mode corresponding to the phase differ- 
ence between two condensates, it has been predicted that 
it also binds fractional vortices winding around two dif- 
ferent components by a sine-Gordon kink [6j, resulting 
in a two- vortex molecule, a dimer [9j. However, such a 
molecule structure has not yet been observed in exotic su- 
perconductors except for indirect evidences [32]. In two- 
component BECs of atoms with two hyperfine states such 
as 87 Rb [28J, two condensates can be coherently coupled 
by introducing a Rabi oscillation, which gives the same 
interaction as the Josephson coupling in superconductors. 
In fact, with the Rabi couplings, vortices winding around 
two (or more) different condensates are found to consti- 
tute a dimer [15j[T8] (or a trimer [21]), connected by one 
(or more) sine-Gordon kink(s) [12J. The molecule struc- 
tures are more accessible in BECs than superconductors 
because of stronger repulsion between two kinds of vor- 
tices and tunable atomic interactions in experiments in 
BECs. However, effects on this term were not studied in 
the vortex phase diagram of two-component BECs. 

The purpose of this Letter is to systematically study 
effects of the Rabi coupling in vortex lattices in two- 
component BECs under rotation. In the limit of a strong 
Rabi coupling, each vortex molecule is tightly bound to 
become an integer vortex, where one can expect the usual 
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Abrikosov lattice of integer vortices. The effects of the 
coherent coupling were simulated in [33J for an atomic- 
molecular BECs mixture. However, the richer vortex lat- 
tice structure in atomic two-component BEC gives rise 
to a wider variety of configurations, as shown below. A 
quite nontrivial question is how the vortex lattices with- 
out the Rabi coupling known before are deformed into 
the integer Abrikosov lattice when the Rabi coupling is 
gradually increased. We find various new structures of 
vortex lattices and new phenomena of vortices, i.e., a 
bound state of multiple dimers bound by an intermolec- 
ular force and exchanging partners among multiple multi- 
dimer bound states analogous to chemical reactions. For 
small Rabi couplings, fractional vortices in triangular or 
square lattice constitute vortex hexamers or tetramers, 
respectively. With increasing the Rabi coupling, they are 
broken into a set of dimers. Then, they exchange their 
partners to organize themselves to prepare for becoming 
the Abrikosov lattice of integer vortices. 

The energy functional of the Gross-Pitaevskii equa- 
tions for the rotating BECs subject to a trapping har- 
monic potential can be written as: 
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where the derivatives include (r, 0) coordinates. We mea- 
sure distances and energies in terms of b = ^/h/muo and 
huj respectively, where m is the mass of the atoms and u 
is the frequency of the trapping harmonic potential. 

The phase diagram of the vortex lattice forming in the 
condensate was studied in [131 [14] and a rich variety of 
lattices was found. The structure of the two component 
lattice depends on the sign and on the magnitude of #12, 
the coupling constant of the intercomponent interactions. 

When #i2 /# = S < 0, the vortices of different compo- 
nents are attracted and are combined into integer vor- 
tices. Because of repulsion among them, they organize 
in a triangular lattice. If S > 0, depending on the value 
of S and vortices in each component organize in the 
triangular lattice, the square lattice, or the vortex-sheet. 
When 5 = the vortices of each component are organized 
in an Abrikosov lattice, but the vortices of different com- 
ponent are decoupled. If 5 is increased, the intercom- 
ponent interaction results in a repulsive force between 
the vortices of different components [THl HO] and the two 
component lattice has an hexagonal structure. When 5 
is increased further the unitary cell of the single compo- 
nent lattice is changed from a triangle to a square. The 
value of 5 for which the lattice reorganize depends on 
rotation speed. If 5 > 1, phase separation occurs and 



1/10 




Integer vortex Abrikosov lattice 



t t 1 



I 3(a) || 3(b) || 3(c) || 2(g) 

i/4 t X 



2(e) 2(f) 



2(h) 



2(a) 2(b) 2(c) 2(d) 



ao 




bo 




CO 




do 



5 = 912/9 



0.3 



FIG. 1: The vortex lattice phase diagram, ao, 60, Co, do indi- 
cate the lattices found in [14 ]. ao corresponds to the triangular 
lattice, bo to the square lattice, and do to the vortex sheet, 
while Co represents the configuration staying at the boundary 
between lattice and vortex sheet regions. The others are the 
different configurations explained in the text. The labels refer 
to the pictures of Fig. 2 and Fig. 3. 



vortex sheets appear. 

We now want to systematically study the effect of an 
internal coherent coupling on the vortex lattice. When 
the two different components are atoms with different hy- 
perfine spin states, the coherent coupling can be achieved 
by Rabi oscillations. This results in the formation of a 
two vortex molecule, namely a dimer. 

The potential induced by Rabi oscillations has the form 



Vk = -w*(***2 + *2*i). 



(2) 



This interaction couples the atoms of the two components 
via the relative phase of the order parameters. We expect 
the interplay between the Rabi interaction and the atom- 
atom interaction triggered by #12 to be fundamental for 
vortex lattice structure. Then we consider the ratio uur/5 
as the relevant quantity. 

We first reproduced the results of [14J by numerically 
minimizing the free energy 0. They are reported on the 
horizontal axis of Fig.[l] We minimize the free energy 
by the non-linear conjugate gradient method (the imagi- 
nary time propagation) in the FreeFem++ package. We 
calculate the ground state of the system with ujr > by 
minimizing the energy functional ([I]) with ([2| added. We 
use the converged lattice solution obtained with ujr = 
as the starting-point configuration and then we increase 
the Rabi frequency by steps Aur = 5/20. We then take 
the results of the converged calculation obtained for some 
value of ujr as initial configuration for the next step. 
By using this "adiabatic" evolution method we are confi- 
dent that we always obtain the true ground state for the 
vortex lattice. In all the numerical simulations we take 
9i = 92 = 9 = 1, & = 0.98 and /ii = /i2 = 15. The di- 
agram obtained is schematically reported in Fig. [T] For 
simplicity we reported labels for the various configura- 
tions which we explain below. 



3 



(a) S 

15 



0.2,^ 0.03 



in 



(b) (5 = 0.5, a;* = 0.075 

15 r 

10 



A A * • • 
***** 

• • • •/ 




-15 



10 15 



-15- 



(e) 8 

15 



10 -5 5 10 15 

= 0.2, u)r = 0.05 




-15-10-5 



10 15 



(f) 5 = 0.b,UR = 0.125 

15 r 



-15-10 -5 5 10 15 





G-O G-O 

G-O G-O 



-10 -5 5 10 15 



-15-10 -5 5 10 15 



(c)S = 

15 



l,u R = 0.05 (d) 6 = 1.1, ujr = 0.11 

15 r 




-15-10 -5 5 10 15 -15-10 -5 5 10 15 



(g) 6 = 1,ur = 0.25 (h) 5 = l.l,u R = 0.275 
15 i 1 15 



' • • • «» l 



^ _l ^^^^^ 



-15-10 




-15-10 -5 5 10 15 -15-10 -5 5 10 15 



FIG. 2: (color online) In all subfigures the top panel is a plot of the Rabi energy in Eq. Q (white is positive, identified with 
domain walls) and the bottom panel is a plot of the density profile of the condensate n— |^i| 2 + |^2| 2 (black dots are vortices 
in either of the two components). The inset is a schematic drawing of the lattice structure, where the two components are 
distinguished by empty or filled circles, which are joined by domain walls depicted by a red line. The values of the parameters 
5 and ujr are shown for each case. In (a) a lattice of hexamers is visible, while in (b) we find tetramers. The tetramers are 
slanted in (c) with a different value of 5. No multi-bound states of dimers appear for 5 > 1. The lattice defect in (d) appears 
just by chance and in fact it resolves with increasing ujr, as can be seen in subfigure (h). Panels (e-g) show the same cases of 
(a-c) with higher ujr. In these cases the bound states are broken into dimers. To visualize the domain walls, the plots of the 
phase difference between the two condensates are given in Fig. [4] in the appendix. 



We find that for ujr/S < 1/10 multi-dimer bound 
states appear. When 6 takes values corresponding to 
triangular lattices < 6 < 0.3, the system exhibits 
vortex hexamers connected by domain walls, as shown 
in Fig. [2|a); if 5 is within the square lattice range 
0-3 ^ 8 < 1, each bound state is composed of four vor- 
tices making up tetramers, as shown in Fig. 2(b); when 
5 = 1 slanted tetramers appear, as shown in Fig. 2(c). 
This behavior is signaled by the form of the Rabi energy, 
whose maxima lie between a group of vortices. When the 
tension of neighboring domain wall and anti-domain wall 
is small enough they can bend toward each other and con- 
stitute a bound state, because of the attraction between 



them. The structure of the lattice is schematically re- 
ported in the inset in the bottom panel of Figs.|2](a),|2|b) 
and|2jc). However, these bound states are metastable 
and when 1/10 < ljr/S < 1/4 they split into dimers, as 
in Figs. [2|e), |2|f), and 2(g). The phase separation for 
5 > 1 prevents from the formation of multi-dimer bound 
states and only single dimers are formed even for higher 
values of ujr/5, as shown in Fig. [2^d) and [2jh). This 
happens because, for S > 1, even with very small ujr, 
there are no dimers whose domain walls are so close to 
bend to form multi-dimer molecules. This can be seen 
from Fig.^d) and^h). 

When ujr increases further, the lattice modifies dras- 



4 



6 = 0.2, uj R = 0.06 



T4L 



-15-10 -5 



5 10 15 




-15-10-5 



10 15 



(b) 5 = 0.5, u;* = 0.15 
15 



(c) S = 0.9, ujr = 0, 
15 



315 




;g^6':(^):c 



OG 3© OB 

8 % OO O-O OrO ee e© e© e© e 

sa g«> CQ~^XQ~Oa Q~~Q) ©© ©© o© q© e 

gs f« O*0"^0""<^€> 66 99 00 09 69 

®® e© e© ©e e© e 

H (^OX§^)(Q^G) o© o© s© e© q© 

n ™ o€> ^GToTeTJ 



10 15 



g© e© ©e 



e© ©© o© 
3 e© g© e© 
oe e© e© 
S oo oc 
o© a© e© 
s e© o© e© 
o© ©© a© 



FIG. 3: (color online) The different partner changing patterns that we have identified. The top panel in each subfigure shows 
the plot of the Rabi energy in Eq. |2) (white is positive, identified as domain walls) and the middle panel shows the density 
profile of the condensate n— |^i| 2 + ]^2| 2 (black dots are vortices in either of the two components). Schematic reproductions 
of the modifications on the dimers and of the lattice itself after the partner changing took place are shown at the bottom of 
each column. The two components are distinguished by empty or filled circles, joined by domain walls depicted by a solid 
straight line. In (a) we show the partner changing process for the hexagonal lattice, appearing when S < 0.3. The domain walls 
binding the vortices are broken and reconnected to be new sets of dimers indicated by the dashed round boxes; the dimers 
indicated by solid round boxes are instead kept unchanged. The defect in (a) is physically relevant, representing a domain wall 
junction between different domains. The subfigures (b) and (c) represent the two partner- changing patterns realized in square 
vortex lattice, 0.3 < S < 1, with the same conventions of picture (a). In these two cases no defect appeared in the lattice, which 
exhibits a periodic structure. To visualize the domain walls, the plots of the phase difference between the two condensates are 
given Fig. [5] in the appendix. 



tically. The reason for this modification is that for large 
values of ujr the dimers are undistinguishable from inte- 
ger vortices, in which the positions of the vortices of dif- 
ferent components coincide. Integer vortices repel each 
other and are organized in an Abrikosov lattice. How- 
ever to reach the Abrikosov configuration, the various 
lattices of fractional vortices must be deeply modified. 
This modification is achieved when 1/4 < uor/5 < 1/3. 
The vortices change their partners. 

We were able to find the elementary patterns of re- 
arrangement, which we depicted in Fig. [3j In the bot- 
tom row of each subfigure schematic representations of 
the process of partner changing are reported. We indi- 
cate with round solid boxes the dimers left unchanged, 
while the dashed boxes represent the formation of new 
dimers. When < S < 0.3 the lattice is triangular for 
< 1/4, but when uor is increased the domain walls 
connecting vortices in dimers can break and reconnect as 
depicted in Fig. |3ja). The partner changing process can 
be realized along three different directions. This is due to 



the discrete rotational symmetry of the two-component 
hexagonal vortex lattice appearing when 5 < 0.3. In the 
calculation reported in Fig. [3^a) all these possibilities are 
realized, dividing the lattice in three domains separated 
by a domain wall junction, as can be seen from the pic- 
ture in the middle panel where the density profile of the 
condensate is shown. 

If 0.3 < 5 < 1, the vortices can change partner in two 
different ways, as shown in Fig. [3^b) and|3jc). The bot- 
tom panels represent the patterns of the partner chang- 
ing processes, with the same conventions of Fig.ffla). As 
can be seen from the schematic diagram in Fig.lSFb), four 
vortices forming a square cell change their partner, while 
links in the adjacent square cell are kept unchanged. The 
situation is different in Fig. [3^c), where vortices change 
their partner in alternating rows. The disposition of the 
dimers depicted in the bottom right panel of Fig. [3^c) is 
similar to that appearing in Fig. [2|h). However the rel- 
ative ordering of the vortices belonging to different com- 
ponents is different. This can be understood by recalling 
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that for 5 > 1 the free energy is minimized when vortices 
of the same component form vortex sheets. Instead in 
the case of Fig. |3|c), the lower energy is achieved when 
vortices belonging to the same component are well sepa- 
rated. 

In conclusion we have found a rich variety of lattices 
which can appear in two-component BECs where internal 
coherent coupling is induced. We introduced the inter- 
action via Rabi oscillations and studied systematically 
its effect on the vortex lattice, increasing gradually the 
frequency uor from uor = 0. The vortex lattice is reorga- 
nized when ur is increased up to high values in order to 
reach the triangular lattice configuration where all dimers 
become integer vortices. During this process multi-dimer 
bound states are formed and vortices can exchange part- 
ners in different ways, depending on the form of the frac- 
tional vortex lattice. The patterns of the exotic vortex 
lattices which we have found in this Letter can be easily 
distinguished from usual triangular or square lattice from 
observations for instance by the time of flight. 
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Appendix 

Here we present more results about the domain walls 
connecting the vortices when coherent coupling is in- 
duced between two components. We want to show more 
specifically the bending which the domain walls exhibit 
when the Rabi coupling ujr is small and multi-dimer 
bound states are formed and also clarify more the part- 
ner changing schemes. If we take the wave functions of 
the two component with the form 

*x(r) = |* i(r)|e iei W , * 2 (r) = \H> 2 (r)\e i0 ^ , (3) 

we can define the phase difference between the two com- 
ponents as 

<P(r)=e 1 (r)-0 2 (r). (4) 

We can then write the Rabi energy term as 

V R = -o; fl (*I(r)* 2 (r) + tf£(r)tf i(r)) 

= -2o; fl |*i(r)||* 2 (r)|cos(0(r)). (5) 

To decrease this term, 0(r) must be close to zero mod- 
ulo 27r, with the amplitudes of the wave functions in- 
creased in the overlapping region. Then sine-Gordon do- 
main walls form between vortices. 

In our calculation we see that this domain wall can 
bend and then multi-dimer bound states can form in 
the vortex lattice, when the Rabi coupling ujr is small 
enough. In Fig. [4] we show the plots of the phase differ- 
ence (j)(r) for some set of the parameters {5 = gi2/g, &r}. 
Fig. |4ja)-(h) correspond to FIG. 2 (a)-(h) in the main 
text, respectively. In the top panels we plot ([5| where 
the whiter regions are maxima. The middle panels are 
plots of the phase difference, in which domain walls be- 
tween vortices are clearly visible. In the bottom panels 



we superimpose the phase difference plot with the plot 
of the condensate density n = |^i| 2 + 1^2 1 2 , where the 
light gray dots at the extrema of domain walls are vor- 
tices belonging to different components. In the first three 
columns we see that domain walls are bent and multi- 
dimer molecules are formed. In the case of Fig. |4[a) we 
find hexamers, while tetramers appear in Fig. RTb). In 
Fig. [4^c) domain walls are less bent than in the two pre- 
vious cases because the elementary cell of the lattice is 
a slanted square, so that it is more difficult for them to 
interact with each other. In all the three cases, we notice 
from the Rabi energy plot that the domain walls inter- 
act and bind together three or two dimers. In Fig. |4^d) 
domain walls are not bent at all, even with very small 
Rabi coupling. This is due to the fact that the dimers 
are disposed in an alternating pattern, so that no domain 
wall is faced with the other. 

In Fig. [4] (e)-(h) we report the same plots obtained 
with the same values for 5, but with higher ujr. We can 
see that the domain walls are now mostly straight and 
all the bound states split into single dimers. 

In Fig. [5] we show the analogous plots for a higher 
particular value for ujr, at which vortices change their 
partner with which they form the dimers, as explained 
in the main text. Fig. [5ja)-(c) correspond to FIG. 3 
(a)-(c) in the main text, respectively. From this plots 
it is clear the disposition of the dimers after the partner 
changing process. In Fig. [5|a) we see that some domain 
walls are bent and multi-dimer bound states are possible 
in the center of the cloud, where a domain wall junction 
is created between the three different domains in which 
the lattice is divided. Fig. [5jb) shows that the dimers 
are now formed with different orientation, compared with 
Fig.Qf); Fig.[5|c) presents a different pattern of partner 
changing, realized when S is close to one. 
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FIG. 4: In all subfigures the top panel is a plot of the Rabi energy in Eq. |5| (white is positive), the middle panel is a plot of 
the phase difference <f>(r) where the domain walls connecting vortices are clearly visible and the bottom panel is a superposition 
of the middle panel with the plot of the condensate density n — |^i| 2 + |^2| 2 (light gray dots are vortices in either of the two 
components). The values of the parameters 5 and ujr are shown for each case. In (a) a lattice of hexamers is visible, while in 
(b) we find tetramers. The tetramers are slanted in (c) with a different value of 5. No multi-bound states of dimers appears 
for S > 1, in (d). In (a), (b) and (c) domain walls are bent and multi-dimers bound states are formed. In (e), (f) and (g) the 
bound states are mostly split into single dimers. In (h) we see that the lattice is periodic with straight domain walls. 
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(a) 5 = 0.2, u R = 0.06 (b) S = 0.5, u R = 0.15 (c) S = 0.9, u R = 0.36 
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FIG. 5: The different partner changing patterns that we have identified. In all subfigures, the top panel is a plot of the Rabi 
energy in Eq. |5]) (white is positive), the middle panel is a plot of the phase difference 0(r) where the domain walls connecting 
vortices are clearly visible, and the bottom panel is a superposition of the middle panel with the plot of the condensate density 
n — |^i | 2 + 1 2 (light gray dots are vortices in either of the two components). In (a) we show the partner changing process 
for the hexagonal lattice, appearing when 5 < 0.3. The defect in (a) is physically relevant, representing a junction between 
different domains. The subfigures (b) and (c) represent the two partner-changing patterns realized in square vortex lattice, 
0.3 < 5 < 1. In these two cases no defect appeared in the lattice, which exhibits a periodic structure. 



